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Primordial black holes (PBHs) are gravitationally collapsed objects 
that may have been created by density fluctuations in the early universe 
and could have arbitrarily small masses down to the Planck scale. Hawk¬ 
ing showed that due to quantum effects, a black hole has a temperature 
inversely proportional to its mass and will emit all energetically allowed 
species of fundamental particles thermally. PBHs with initial masses of 
order 5x10^^ g should be expiring in the present epoch with bursts of 
high-energy particles, including gamma radiation in the GeV-TeV energy 
range. The Milagro high-energy observatory, which operated from 2000 
to 2008, is sensitive to the high end of the PBH evaporation gamma-ray 
spectrum. A search of five years of Milagro data set a local (parsec-scale) 
upper limit of 3.6x10"^ PBH bursts/year/pc3. We will also report the sen¬ 
sitivity of the High-Altitude Water-Cherenkov (HAWC) observatory to 
PBH evaporation events. Finally, we investigate the final few seconds of 
black hole evaporation using Standard Model physics and calculate energy 
dependent PBH burst time profiles in GeV/TeV range. We calculate PBH 
burst light curves observable by HAWC and explore search methods and 
potentially unique observational signatures of PBH bursts. 
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1 Introduction 


Primordial black holes (PBHs) are black holes created in the early universe. Depend¬ 
ing on the formation scenario, the PBHs could have masses ranging from the Planck 
mass to more than million solar masses [1]. In 1974, Hawking showed by convolving 
general relativity, thermodynamics and quantum field theory that a Black Hole (BH) 
has a temperature inversely proportional to its mass and radiates with thermal spectra 
photons and massive particles [2] . As the BH emits this radiation, its mass decreases 
and hence its temperature and emission rate increase. A PBH that formed with an 
initial mass of ~ 5.0 x 10^^ kg in the early universe should be expiring today [3] with 
a burst of high-energy particles, including gamma-rays in the MeV to TeV energy 
range. Thus PBHs are candidate gamma-ray burst (GRB) progenitors [1]. 

Conhrmed detection of a PBH evaporation event would provide valuable insights 
into many areas of physics including the early universe, high energy particle physics 
and the convolution of gravitation with thermodynamics. Conversely, non-detection 
of PBH evaporation events in sky searches would place important limits on models 
of the early universe. 

The properties of the PBH final burst depend on the physics governing the produc¬ 
tion and decay of high-energy particles. In the Standard Evaporation Model (SEM) 
which incorporates the Standard Model of particle physics, a BH should directly 
Hawking-radiate the fundamental Standard Model particles whose Compton wave¬ 
lengths are of the order of the black hole size [5]. As the BH evaporates and loses 
mass over its lifetime, its temperature surpasses the rest mass thresholds of further 
fundamental particle species. When the temperature exceeds the Quantum Chromo¬ 
dynamics (QCD) confinement scale (~ 300 MeV), the BH should directly emit quarks 
and gluons which fragment and hadronize (analogous to jets seen in high-energy col¬ 
lisions in terrestrial accelerators) as they stream away from the BH into the particles 
which are stable on astrophysical timescales EE]. Thus according to the SEM, the 
evaporating black hole will be seen astronomically as a burst of photons, neutrinos, 
electrons, positrons, protons and anti-protons. 

2 Photons from a PBH Burst 

2.1 Hawking Radiation 

Hawking showed that a black hole radiates fundamental particle species at an emission 
rate of O |T] 

(PN _ r/ 27 rh 
dM ~ e- - (-1)2- 

where s is the particle spin and Udof is the number of degrees of freedom of the particle 
species (e.g. spin, electric charge, color and flavor). The dimensionless quantity x is 
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defined by 


E 


( 2 ) 


_ SttGME 

hc^ kTsH 

where E is the energy of the emitted particle, M is the black hole mass, Tbh oc 1/M 
is the black hole temperature, and G, c, and h are the universal gravitational con¬ 
stant, speed of light, and the reduced Planck constant, respectively. The absorption 
coefficient P depends on M, and s. For an emitted species of rest mass m, P at 
E 3 > mc^ has the form 

r(M.B,s) = 27(^)k,W (3) 

such that 7 s (x) —)■ 1 for large x. 

To calculate the spectrum of the hnal photon burst from the PBH, two important 
relations pertaining to the hnal phase of BH evaporation are needed. The hrst relation 
is the black hole mass M expressed as a function of remaining lifetime. The mass 
loss rate can be written as |1] 

dM _ a{M) 

where the factor a{M) incorporates all emitted particle species and degrees of free¬ 
dom. As the BH evaporates, the value of M is reduced by an amount equal to the 
total mass-energy of the emitted particles. By conservation of energy. 


dt ^ Jq dEdt 


(5) 


where the summation over i is over all radiated species. Therefore, 


a(M) 



dEdt 


EdE. 


( 6 ) 


Assuming the Standard Model including the top quark, 125 GeV Higgs, and 3 families 
of massive Majorana neutrinos, we hnd that the asymptotic value osm of a{M) as M 
decreases is 

asM = 8.40 X lO^’^kg^s"^ (7) 

For the current and future generations of high energy gamma-ray observatories, we 
are interested in bursts generated by black holes of temperature Tbh ^ ^ TeV. Thus 
for Tbh ^ 1 TeV BHs (corresponding to M < 10^ kg and a remaining evaporation 
lifetime of r < 500s), a{M) ^ asM and the BH mass as a function of remaining 
evaporation lifetime r is 

M(r) (3asMn) ^ = 1.36 X 10® kg. (8) 
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The second relation that we require for the hnal evaporation phase is the temper¬ 
ature Tbh expressed as a function of r. Combining Equations and we have 

for kTsH > 1 TeV. 

2.2 QCD Fragmentation 

According to the SEM, Equation [T] describes the direct Hawking radiation of the fun¬ 
damental particle species of the Standard Model: the leptons, quarks, and the gauge 
bosons As they stream away from the BH, these particles will then evolve by 

Standard Model processes, ultimately into the particles which are stable on astrophys- 
ical timescales: photons, neutrinos, electrons, positrons, protons, and antiprotons. In 
particular, quarks and gluons will undergo fragmentation and hadronization into in¬ 
termediate states which will eventually decay into the astrophysically stable particles. 

For application to PBH searches at gamma-ray observatories, we seek the total 
photon emission rate from the BH. The photon spectrum has several components: (i) 
the “direct photons” which are directly Hawking radiated by the BH: this component 
peaks at a few times Tbh and is most important at the highest photon energies at 
any given Tbh', (h) the “fragmentation photons” arising from the fragmentation and 
hadronization of the directly Hawking radiated quarks and gluons: this component 
is the dominant source of photons at energies below Tbh', and (hi) the photons pro¬ 
duced by the decays of other Hawking-radiated fundamental particles, e.g. the tan 
lepton, W and Z gauge bosons, and Higgs bosons: this component is small compared 
to the component produced by the directly Hawking radiated quarks and gluons and 
is neglected here. (We note that, because the W, Z, and Higgs bosons decay predom¬ 
inantly via hadronic channels, the main effect of component (iii) will be to enhance 
the fragmentation photon component (ii) somewhat.) 

2.3 The Pion Fragmentation Model 

In jet fragmentation and hadronization, most of the photons arise from the decays of 
7r° states, which is decays to 27 ’s with a branching fraction 98.8%. We proceed 
assuming that the QCD fragmentation of quarks and gluons may be approximated 
entirely by pion production. In our pion fragmentation model, two questions must be 
addressed: what is the pion distribution generated by the partons (the initial quarks 
and gluons) and what is the photon spectrum generated by the pion decays? 

To answer the hrst question, we utilize a heuristic fragmentation function 

(1 - zf ( 10 ) 


3 



where z = Ej^/E is the energy fraction carried by a pion generated by a parton of 
energy E |1]. We assume the same form of DT^ji^z) for all partons and that all of 
the parton energy goes into pions i.e. the function in Equation is normalized such 
that Jg zDT^/i{z)dz = 1. The instantaneous BH pion production rate per pion energy 
interval is then 


dE-i^dt 



d^N, 

dEdt 




zE)dzdE. 


( 11 ) 


2.4 Photon Flux from Pion Fragmentation 

To answer the second question, we obtain the photon flux from the —)■ 27 decay of 
the pion distribution. Because the fragmentation function includes all three 

pion charge states 7 r''', 7 r“, and 7 r° as equal components, and each 7 r° decays into two 
photons, we must multiply by 2/3 to get the 7 multiplicity. In the 7 r° rest frame, 
the two photons have equal energies, m^c^/2 and equal but opposite momenta. In 
the reference frame of the gamma-ray detector, the energies of the two photons are 
unequal but complementary fractions of the vr^ energy in the detector frame, We 
assume that only one of the photons in each pair is detected]^ 

Let 9' be the angle between the momentum of the observed photon in the 7 r° rest 
frame and the momentum of the pion in the detector frame. In the detector frame, 
the photon energy E-y = {E.„/2){1 + (3 cos 9') where the 7 r° velocity (3 = v/c ~ 1 and 
Et, = tUtt/ y/l — [3“^. Because the angular distribution of the photons is isotropic in 
the 7 r° rest frame, the distribution of photon energy in the detector frame is 


d^N^ 

dE^dt 


2 2'Kdcos9' /•“ 

3 y_i dvr dE^dt ^ 


{E^/2){1 + (3 cos 9')]dE^. 


( 12 ) 


2.5 Parametrization of the Direct and the Fragmentation 
Contribution 

To further simplify our analysis and subsequent calculations, we parameterize the di¬ 
rectly Hawking radiated and the pion fragmentation components of the instantaneous 
PBH burst photon spectrum using the variable 

5,.1,287x10-‘(^)(H)‘'“ (13) 


as follows: 

* The angle between the 2 photon trajectories in the detector frame will be very small because 
of the large Lorentz boost. However, if the BH is at a distance of order 1 parsec from the detector, 
then only one of the photons from each 7r° decay will hit the detector. 
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f = 1-13 X 10^®(^^)®(exp(^^) - 1) ^ X F{Q GeV ^ s ^ (14) 

\dh^dtj 

where 

^(^7) = 1-0 ^7 < 2 

F{i^) = exp([-0.0962 - 1.982(ln^^ - 1.908)][1 + tanh[20(ln^^ - 1.908)]]) for > 2 
and 

(;S^) " 6.339 X 10"p?7-’''"|l - BsKt - 0.3)] 

V '^'^7 ^ / frag. 

+ 1.1367 X 10^4 exp(-e.^)[^.^(^.^ + 1)]"^ Qs{i^ - 0.3) GeV^ st*) 


where 


©^(A) = 0.5[1 + tanh(10 A)]. 


These parametrizations are valid for energies above ^ 1 GeV and are accnrate to 
± 15 % for in the range from 0.1 to 5.0, and to ± 3 % for smaller and larger 
If greater accnracy is reqnired, we nse a table of the ratios of the exact integration 
of Eqnation IT using the fragmentation function 1^ to the parameterized value, to 
correct the parameterized values. 


2.6 PBH Burst Light Curve 

To hnd the time evolution of the PBH burst, we integrate the total instantaneous 
photon flux d^N ^/ (dE-ydt) over photon energy Ey while retaining the time depen¬ 
dence: 


dNy 

dt 


Emission 



4^dEy. 

dEydt 


(16) 


In general Emin and Emax are set by the energy range of the detector. As a case study, 
we use here the High Altitude Water Gherenkov (HAWG) observatory [B] for which 
the relevant energy range is Emin=50 GeV to E max =100 TeV. 

The time prohle (the light curve) of the PBH burst at the detector can be calcu¬ 
lated as follows: 


dNy 

dt 


Detection 


dvrr^ 


T(E 




dE^ 


^^dEydt ^ 


(17) 
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Figure 1: Simulated PBH Burst lightcurve observed by HAWC (at a distance of 0.05 
light years) obtained by convolving with the HAWC effective area published in Ref [7] . 
This shape is also well described by a power law with a index of -0.66. 



Figure 2: Simulated detection time prohles of a PBH burst at a distance of 0.05 light 
years observed by HAWC in multiple energy bands. Photon numbers detected in each 
energy band is shown in the legend. 
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where A{E^) is the effective area of the HAWC detector and r is the distance to the 
PBH. Figure shows the time profile of a PBH burst at a distance of 0.05 light years 
in the HAWC energy range and for the HAWC effective area published in Ref [7] . 

It is also interesting to investigate the energy dependence of the PBH burst time 
prohle. In Figure]^ we show dN-y/dt calculated using various energy bands in the 
HAWC energy range. As seen in Figure the low energy bands show similar emission 
prohles. However, above ~ 10 TeV the burst emission time profile is energy- 
dependent and exhibits an inflection around r ~ 0.1 seconds. This occurs due to the 
domination of the direct photon component at the highest energies. 

3 Summary and Conclusions 

In this paper, we have briefly reviewed the theoretical framework of the PBH standard 
emission model and calculated the PBH burst light curves which would be observed 
by HAWC. The main hndings and conclusions of our paper are: 

1. We have developed approximate analytical formulae for the instantaneous PBH 
spectrum which includes both the directly Hawking radiated photons and the 
photons arising from the other directly Hawking radiated species. 

2. For the first time, we have calculated the PBH burst light curve and studied its 
energy dependence at a detector. 

3. For low energies {E^ < 10 TeV) the light curve prohle does not show much 

evolution with energy and is well described by a power law However, at 

high energy the light curve displays signihcant energy dependence that may be 
used as an unique signature of PBH bursts. The HAWC observatory is sensitive 
in this high energy range and potentially can be used to uniquely identify PBH 
bursts. 
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